Experimental filtering of two-, four-, and six-photon singlets from single PDC source 
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Invariant entangled states remain unchanged under simultaneous identical unitary transformations 
of all their subsystems. We experimentally generate and characterize such invariant two-, four-, and 
six-photon polarization entangled states. This is done only with a suitable filtering procedure of 
multiple emissions of entangled photon pairs from a single source, without any interferometric over- 
laps. We get the desired states utilizing bosonic emission enhancement due to indistinguishability. 
The setup is very stable, and gives high interference contrasts. Thus, the process is a very likely 
candidate for various photonic demonstrations of quantum information protocols. 
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Entanglement is an essential tool in many quantum 
information tasks. Entangled states of two qubits proved 
to be useful in various quantum communication protocols 
like quantum teleportation, quantum dense coding, and 
quantum cryptography. They are the essence of the first 
versions of Bell's theorem[l|. However, the expansion of 
quantum information science has now reached a state in 
which many schemes are involving multiparty processes, 
and could require multiqubit entanglement. 

There is an interesting series of multiqubit states, 
l^jt ), where k = 2,4,6 or more. They are invariant 
under actions consisting of identical unitary transfor- 
mations of each qubit [3]: C/® fc |^~) = ), where 
U® k — U ® ... ® U denotes a tensor product of k iden- 
tical unitary operators U. The property protects the 
states against collective noise. The states are useful e.g. 
for communication of quantum information between ob- 
servers who do not share a common reference frame 
any realignment of the receiver's reference frame corre- 
sponds to an application of the same transformation to 
each of the sent qubits. The states | Vf^ ) can also be used 
for secure quantum multiparty cryptographic protocols, 
such as the multi-party secret sharing protocol ToL 11 1. 

We generate correlations which characterize the six- 
photon | ) entangled state. This is done in a six pho- 
ton interference experiment. A six-photon interference 
was reported recently in [13] • To obtain graphs states the 
authors of ref. 12:] used three pulse-pumped paramet- 
ric down-conversion (PDC) crystals, and interferometric 
overlaps, to entangle independently emitted pairs (each 
from a different crystal) with each other. Schemes of this 
kind are generalizations of those of ref. [13]. However, 
the overlaps make the scheme fragile. 

In our experiment, by pulse pumping just one crystal 
and extracting the right order process via suitable filter- 
ing and beamsplitting (the method of [IH), we observe 
simultaneously effects attributable to the multi-photon 
invariant entangled states | ^ ) , | J ) , and | ^ ) . The 
setup has no overlaps and therefore no interferometric 



alignment is needed. It is strongly robust, and the out- 
put is of high fidelity with respect to the theoretical states 

I**)-. 

A simple quantum optical description of two phase 
matched modes, of the multiphoton state that results 
out of a single pulse acting on a type-II PDC crystal, can 
be put as 



0V u 0H 



)]"|0). 



(1) 



The symbol a) 0H (b^ v ) denotes a creation operator for one 
horizontal, H, (vertical, V") photon in mode ao (bo), etc. 
C is a normalization constant, the coupling parameter a 
is a function of pump power, non-linearity and length of 
the crystal. This is a good approximation of the actual 
state, provided one collects the photons under conditions 
that allow full indistinguishability between separate two- 
photon emissions [3]. First, second, and third order 
terms in the expansion in eq. |T]), correspond to an emis- 
sion of two, four, and six photons, respectively, into two 
spatial modes. These terms can be re-interpreted as the 
following superpositions of photon number states: 

\lH ao ,lV bo )-\lV ao ,lH bo ), (2) 

1 2H ao , 2V bo ) - | lH ao , IKo , lV bo ,lH bo ) 

+ \2V ao ,2H bo ), (3) 



1 3-ffao ! 3Vb ) - 1 2H ao , lV ao , 2V bo , lH bo ) + 

+ | lH ao , 2V ao , 1 Ho , 2ff 6o ) - 1 3V ao , 3H bo ) , (4) 

where e.g. 2H ao and 3H ao denotes two and three horizon- 
tally polarized photons in mode ao, respectively, etc. The 
second and third order PDC is intrinsically different than 
simple products of two and three entangled pairs. Due 
to the bosonic nature of photons, emissions of completely 
indistinguishable photons are more likely, than the ones 
giving birth to photons with orthogonal polarization. 
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We report a joint observation, in one setup, of the cor- 
relations of the invariant two, four and six-photon polar- 
ization entangled states given by the following superpo- 
sitions: \^) = j=(\HV)-\VH)), 



\^) = ^=\GHZ+)--L\EPR)\EPR), (5) 



and 



\*-) = -L\GHZe) + ^(\W 3 )\W 3 ) -\W 3 )\W 3 )). 



(6) 



The states in the superpositions are given by: 

\GHZ+) = -=QHHVV) + \VVHH))/V2 
V2 



\GHZq) = QHHHVVV] 



-\VVVHHH)), 



and 



\EPR) = -j={\HV) + \VH)). 



Finally 



|W 3 ) = ~^{\HHV) + \HVH) + \VHH)). 
v3 

The ket \ W) is the spin-flipped \ W). The states (??-EJ) 
are obtained out of different orders of the PDC emission 
(figure 1), by selecting specific double, quadruple and 
six-fold coincidences. 

In our setup we use a frequency-doubled Ti:Sapphire 
laser (80 Mhz repetition rate, 140 fs pulse length) yield- 
ing UV pulses with a central wavelength at 390 nm and 
an average power of 1300 mW. The pump beam is fo- 
cused to a 160 /im waist in a 2 mm thick BBO (/3-barium 
borate) crystal. Half wave plates and two 1 mm thick 
BBO crystals are used for compensation of longitudi- 
nal and transversal walk-offs. The third order emission 
of non-collinear type-II PDC is then coupled to single 
mode fibers (SMF), defining the two spatial modes at the 
crossings of the two frequency degenerated PDC emis- 
sion cones. Leaving the fibers the down-conversion light 
passes narrow band (AA = 3 nm) interference filters (F) 
and is split into six spatial modes (a,b,c,d,e,f) by or- 
dinary 50% — 50% beam splitters (BS), followed by bire- 
fringent optics (to compensate phase shifts in the BS's). 
Due to the short pulses, narrow band filters, and sin- 
gle mode fibers the down-converted photons are tempo- 
rally, spectrally, and spatially indistinguishable 14}, see 
Fig. Q] The polarization is being kept by passive fiber 
polarization controllers. Polarization analysis is imple- 
mented by a half wave plate (HWP), a quarter wave plate 
(QWP), and a polarizing beam splitter (PBS) in each 




x PBS 



BS d 



I H > >v 



BS a PBS 1 



FIG. 1: Experimental setup for generating and analyzing the 
six-photon polarization-entangled state. The six photons are 
created in third order PDC processes in a 2 mm thick BBO 
pumped by UV pulses. The intersections of the two cones 
obtained in non-collinear type-II PDC are coupled to single 
mode fibers (SMF) wound in polarization controllers. Nar- 
row band interference filters (F) (AA = 3 nm) serve to re- 
move spectral distinguishability. The coupled spatial modes 
are divided into three modes each by 50%-50% beam split- 
ters (BS). Each mode can be polarization analyzed using half 
wave plates (HWP) and a polarizing beam splitter (PBS). 
Simultaneous detection of six photons (two single photon de- 
tectors for each mode) are being recorded by a twelve channel 
coincidence counter. 



mode. The outputs of the PBS's are lead to single pho- 
ton silicon avalanche photo diodes (APD) through multi 
mode fibers. The APD's electronic responses, following 
photo detections, are being counted by a multi-channel 
coincidence counter with a 3.3 ns time window. The coin- 
cidence counter registers any coincidence event between 
the APD's as well as single detection events. 

The states )(k = 2,4,6) exhibit perfect two, four, 
and six qubit correlations. The correlation function is 
defined as an expectation value of the product of local 
polarization "Pauli" observables. If one limits the mea- 



surement to the local observables cos 9i a 



sin Via, 



(with eigenvectors y/l/2{\L) l ±e ! |i?) ; ) and eigenvalues 
±1), the measurements correspond to linear polarization 
analysis in each spatial mode (I — a, b, c, d, e, /). In such 
a case the quantum prediction for the two photon (in 
modes b and d) correlation function reads: E{9b 1 9d) = 
— cos(0b — 0d)- For the four photon counts (in modes a, 
b, d, and e) correlation function is given by 



E(9 a , 9 b , 9 d , 9 e ) = | cos(# a + 9 b - 9 d - 9 e ) 



+ - cos(9 a - 9 b ) cos(6> c - 9 d ). 



(7) 
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TABLE I: Visibilities of the invariant states | ^ k ) (where k = 
2,4,6) 



k 


Modes 


Visibility 


2 


b,d 


0.962 ± 0.004 


2 


b, e 


0.963 ± 0.006 


2 


o,f 


0.962 ± 0.004 


4 


a, b, d, e 


0.919 ±0.014 


4 


a,b,dj 


0.918 ±0.011 


4 


a,b,e,f 


0.919 ±0.014 


4 


a,b,e,f 


0.919 ±0.014 


4 


b,c,d,f 


0.918 ±0.009 


4 


b,c,e,f 


0.920 ±0.012 


6 


a,b, c, d,e,f 


0.838 ± 0.030 



Finally for the six photon events one has 

E(9 a , 9b, 9c, 9d, 9 e , 9f) — 

1 



i^COs(0 o ±0 b ±0 c ±0 d ±0e±0/), (8) 



where ^ is a sum over all possible sign sequences which 
contain only two positive signs, with the sign sequence in 
the first term, proportional to |, excluded. Due to the 
invariance, the correlation functions for all measurements 
around any single great circle of the Bloch sphere look 
the same. 

Fig.[3]shows three experimentally observed two-photon 
correlation functions, £?(#;,, 6*;), where (I = d,e,f). The 
setting 9b is varied, while the other analyzer is fixed 
at 9i = 9 m = 9 n = it/2. This corresponds to diago- 
nal/antidiagonal, D/A, linear polarization analysis. A 
sinusoidal least-square fit was made to the data. The av- 
erage visibility, defined here, as the average amplitude of 
the three fits, is V 2 = 0.962% ± 0.003%. 

Fig. [3] shows how six experimentally observed four 
photon correlation functions E(9b,9i,9 m ,9 n ) depend on 
9b- The other analyzers were fixed at 9i = 7r/2, where 
(I = a, c), (to = d, e), and (n — e, /). The average value 
of the six visibilities is V 4 = 0.9189% ± 0.0049%. 

Finally, Fig. Q] shows similar data for the ex- 
perimentally observed six photon correlation function 
E(9b,9 a ,9 c ,9d,9 e ,9f). Again 9b was varied with the 
other five analyzers fixed at 9i = tt/2 where (I — 
a, c, d, e, /). The value on the visibility is V — 83.79% ± 
2.98%. 

In table [J we present all the experimentally obtained 
two-, four-, and six- photon visibilities. 

We have compared the observed visibilities with theo- 
retical predictions, Uj. To estimate maximal predictable 
visibilities one can use a less simplified description of the 
two-photon state emitted by SPDC event, and replace in 




FIG. 2: Two-photon polarization correlation function. Modes 
a,d,e and / are analyzed in D/A basis and mode b anal- 
ysis basis is varied around the equator of the Bloch sphere 
(a z cos(8b) + G x sin(6b)). The solid lines show sinusoidal fits 
to the experimental data with a average visibility of V2 = 
0.962% ± 0.003%. 
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FIG. 3: Four-photon polarization correlation functions of 
l^jj). Modes a,c,d,e and / are analyzed in the +/ — basis 
and mode b analysis basis is rotated around the equator of the 
Bloch sphere (a z cos 8b + cr x sin 6b). The figures correspond to 
different implementations of the state, using different modes 
(abde, abdf, abef, bede, bedf and beef). The solid lines show 
sinusoidal fits to the experimental data. The solid lines show 
sinusoidal fits to the experimental data with a average visi- 
bility of V 4 = 0.9189% ± 0.0049%. 
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entangled, if the norm of the N-particle correlation ten- 



FIG. 4: Six-photon polarization correlation function. Modes 
a,c,d,e and / are analyzed in D/A basis and mode b anal- 
ysis basis is varied around the equator of the Bloch sphere 
(azCOs(Bb) + o x sin(9b)). The solid line shows a sinusoidal fit 
to the experimental data with a visibility of 83.79% ± 2.98%. 
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m b lv - a lv b l H ) b y 



eq. 



dt J duj J duJx J dw 2 /(wi)/(w 2 )5(wo)e x 

A(uj q -U3x- uj 2 )(al H (uJi)bl v (uj2) - %v(vi)% H (v2))- 

(9) 

This approach is rich enough to take into account the 
frequency phase matching conditions. The creation op- 
erators depend additionally on frequencies, and obey 
[a x(u),aox>(u')] = Sxx>8xx>(u - w'), etc. The func- 
tion f(u>) represents the shape of the filter transmission 
profiles, and g(u>o) represents the frequency profile of 
the pump pulse. If one chooses f(u>) = exp[—((u>f — 
w)/(2cr/)) 2 ] and g(tu) = exp[— ((uj p — uj) / (2a p )) 2 } where 
cr p , and a/ are the FWHM bandwidth of the pump, and 
the filters, and uif = uj p /2, the following formulas for the 
maximal theoretical visibility as a function of the ratio 
r = <7f/a p can be reached, [l7(. For the four-photon pro- 
cess: y* emp = y/i 4- 2r 2 /(I + r 2 ) and for the six-photon 
process: V^ mp = (l + 2r 2 )/[(l + r 2 /2)(l + 3r 2 /2)]. In our 
experiment we used r exp = AA//(4AA P ) = 0.76. This 
corresponds to V4 = 0.93 and V§ — 0.90. The actual 
measured values of visibility for two, four, and six photon 
interference are very close to the predicted ones, see ta- 
ble |TJ Thus, the fact that our setup use only filtering and 
beamsplitting, has interferometric advantages. In other 
words the obtained four and six particles visibilities are 
almost as high as one get for the ratio r exp . 

The high visibility has the following consequences. 
With just a part of our data one can use the simplest one 
of the "experimentally friendly" entanglement indicators 
introduced in [l8| . It guarantees that N-qubit state is 



sor is higher than 1. For the | <J7 2 
*4 ) again just T x 



^yy ' ^22; for 



we take just T xx , 

Tyyyyi ^zzzz i and. 

for |*g ) the components T xxxxxx , T yyyyyy and T zzzzzz . 
With our data we have obtained the 2.785 ± 0.007, 
2.517 ± 0.011, and 2.29 ± 0.14, respectively for each of 
the case. The entanglement threshold is violated by of 
242, 133, and 9.3 standard deviations. Additionally, ac- 
cording to the criteria given in [p| the state cannot 
be described by a local realistic model, if the sum of 
squares of two out the the listed components is above 
1. This is again achieved by our data: for four parti- 
cles we get T 2 XXX + T 2 yyy = 1.646 ± 0.009 (exceeding 1 
by 74.8 standard deviations) and for six particles we have 
T Lxxxx + Tyyyyyy = 1-52±0.11 (exceeding 1 by 4,5 stan- 
dard deviations). Thus the state can be directly (that is 
without still enhancing its fidelity) utilized in classical 
threshold beating communication complexity protocols 

In summary, we have experimentally demonstrated 
that a suitable filtering procedure applied to a triple 
emission from a single pulsed source of polarization en- 
tangled photons leads to two, four and six-photon, high 
visibility, interference due to entanglement, observable in 
a single setup. We utilize the bosonic emission enhance- 
ment occurring in the emission of three photon-pairs 
in PDC, thus the process is not entirely spontaneous. 
The six qubit state that we observed is invariant with 
respect to simultaneous identical (unitary) transforma- 
tions of all qubits. This makes it particularly useful 
for multiparty quantum communication and general 
quantum computation tasks: it is very robust against 
deformations in transfer. Since we use just one source, 
we avoid alignment problems, and thus the setup is 
very stable. We would like to note that the interference 
contrast is high enough for our source to be used in 
two-, four-, and six party demonstrations of quantum 
reduction of communication complexity in some joint 
computational tasks, and for secret sharing, as well as 
in many other quantum informational technologies. 
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